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Abstract 

By using the adequate modified Priifer variables, precise upper and lower bounds on 
the density of states in the (internal) Lifshitz tails are proven for a ID Anderson model 
with bounded potential. 

There have been numerous rigorous works about Lifshitz tails for the ID-Anderson model 
with bounded potentials (see |KWj for a collection of references). The aim of this note is to give 
a simple proof by passing to the normal form of the transfer matrix at a band edge and then 
using the adequate Priifer variables in this regime. This allows to obtain quite precise estimates 
on the integrated density of states (IDS). The paper concludes with a brief outlook on how and 
why similar techniques lead to perturbative results about the IDS and the Lyapunov exponent 
in between the Lifshitz tails and the center of the band, a regime that will be studied in more 
detail elsewhere. 

1 Result 

Let {H a ) a€ Y, be a family of L-periodic real Jacobi matrices on £ 2 (Z). Each H a is specified by 
2L real numbers t a (l) > and v a (l) where I — 1, . . . , L such that for any state ip £ 2 {Z) 

(H a ip)(n) = -t a (n + l)ijj{n + 1) + v a (n)^(n) - t a {n)i){n - 1) . 

The eigenvalue equation H a ip = Eip for EgI can conveniently be rewritten using the transfer 
matrices G SX(R, 2) over one period: 

f t a (L)^{L) \ E ( t a (l)iP(l) \ 

v ^l-i) ) - ^ \ m ) ■ 

The set S is supposed to be finite here and on it be given a probability measure p = Y2aP^^- 
Of course, p a > and ^2 a Pa = 1. The Tychonov space Q = S xZ x {1, . . . L} is furnished with 
the probability measure P = p® z <S> jYli^i- This measure is of the so-called Bernoulli type 
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and it is invariant and ergodic w.r.t. the natural translation action of Z on Q. Associated 
to each configuration uj = ((er n ) ng z,Z) G Q is a Hamiltonian H w obtained by juxtaposition of 
the periodic blocs according to the configuration uj. More precisely, if n = mL + I + k with 
k — 1, . . . , L, then 

(# w V)(n) = -t ff JHlMn + l)+^(iMn)-t ff JW(n-l): 

with the convention that t Um (L + 1) = t (Tm+1 (l). By construction, (H^wen is a strongly contin- 
uous operator family for which the covariance relation w.r.t. the translations holds |PF| FBelj. 
Note that, even in the trivial case where p is supported by only one point a G S, one has a 
covariant operator family given by the L shifts of the L-periodic operator H a . The spectrum 
of any such covariant family is P-almost surely constant. Using approximate eigenfunctions, 
one easily verifies spec(H u ) C {J a pa ^ Q spec(if CT ). Equality instead of inclusion holds if H a is the 
sum of a periodic background operator and a random potential |PFj . In this situation, a band 
edge of the random operator is a band edge also for one of the periodic operators. 

Two fundemmental objects associated to a covariant family of ID operators are the IDS j\f 
and the Lyapunov exponent |PF| . the focus here being only on the former. One of the equivalent 
definitions of the Stieltjes function j\f is formula (JHJ) in Section El Its support is spec(if^). If 
p = S a , then it is straight-forward to calculate the associated IDS, denoted j\f a , by means of 
Bloch-Floquet theory. 

Let now E v G R be a boundary point of spec(H u ) and spec(H CJ ). Hence E v can be either the 
bottom or the top of the spectrum or the boundary point of an internal spectral gap. Within 
any gap, the IDS j\f(E u ) is constant and therefore it is natural to study 

SAf(e) = ±jV(E u ±e) T M{E V ) , 

where the upper signs are chosen if E v is a lower band edge and the lower signs for an upper 
band edge. As above, one also introduces Sj\f v (e). Close to E v , the IDS Af is very small and its 
scaling is universal and called a Lifshitz tail, honoring the original contribution of L. Pastur's 
teacher. 

Theorem Let E v be a band edge of H v and H^, but not of H a if a ^ v. Suppose moreover 
that the eigenvector of is not an eigenvector ofTf v for any other o ^ v. Then there exist 
constants C > and eo > such that for all e G [0, eo] 

i SK(e) pF™^ +1 < SjV(e) < CSK{e)pF^ . (1) 



Because of the van Hove singularities, this result allows to read off the Lifshitz exponent 

log(|log(&A/"(e))|) = 1 
™ |log(e)| 2 ' 

Moreover it gives precise bounds on the IDS and hence bounds on the so-called Lifshitz constants. 
For such asymptotics to hold, it is crucial that the transfer matrices be uniformly bounded. 
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However, it should be possible to relax the condition that (X, p) be discrete. By approximating 
a continuous measure p by a discrete one, a straight-forward adoption of the presented argument 
allows to obtain at least the lower bound. Finally, the (generic) conditions that E v is only a 
band edge of H v as well as the condition on the eigenvectors of Tf v are not essential, but the 
statement of the result would be a bit more involved. 

As will be discussed in Section EJ the estimates (JQ are in a certain sense based on a determin- 
istic argument, albeit taking place in a random model. The Lifshitz constants characterizing the 
true behavior of the IDS are model-dependent. One approach to calculate them is perturbation 
theory. As will be sketched in Section 0J this also allows to go beyond the Lifshitz tail regime. 

Acknowledgements: This work was supported by the SFB 288. 



2 Setup 

First let us recall some of the analysis of the periodic operators H a . The eigenvalues of the 
transfer matrix Tf are A = \ (Tt{T^) + ^/Tt{T^ ) 2 - 4) and 1/A. Hence, if |Tr(T^)| < 2, there 

are complex conjugate eigenvalues A = e %v and e~ %v and the transfer matrix is conjugate to 
a rotation by the phase r\. This phase is also called the rotation number and one speaks of the 
elliptic case. Then E G spec(H a ) and the IDS J\f a at E is equal to t)/(Ltt) up to a multiple of 
1/L coming from the gap label |.IM[ IBel] . On the other hand, if |Tr(T^)| > 2, the eigenvalues 
are both real. One of them has a modulus bigger than 1 and one smaller than 1. This is the 

hyperbolic case, E ^ spec(if CT ) and the transfer matrix is conjugate to the dilation ( ^ 

Again due to the gap labelling, one has M a {E) = m/L for some positive integer m < L. 

Let now E v be a band edge of the operator H u . One then has |Tr(T^*') | = 2 and A = ±1 and 
is therefore in the so-called parabolic case. The transfer matrix T„ v has only one eigenvector 
denoted v as well as a principal vector w satisfying (T„ v — Xl)w = v. Hence the basis change 

with M = (v w) conjugates T^ v to a Jordan normal form M~ X T^ V M — \ \ \ J • As the 



x A , 

energy varies around E v , one is in either of the above elliptic or hyperbolic cases. However, the 
corresponding basis changes become singular at E v and it is better to work with a basis change 
into an object close to the parabolic normal form. It would be possible to simply work with the 
energy independent M, but for sake of more explicit formulas later on let us choose (which is 
easily seen to be possible for |e| < eo for some eo) an energy dependent basis change M t such 



M -i T E„+e Me / ^ e; ; , n t = 2 - A Tr(T^ +e ) . (2) 



that 

-Ir^Ev+e-Af _ I A(l — « e ) 1 \ 

-K e A 

As T ! f" +e is a polynomial in e, M t is analytic. If n t < 0, one is in the hyperbolic case and 
for n e > in the elliptic one. In the latter the rotation number rj is then given by e lT> = 

+ i^/4:K e — Tcf J . Band touching happens if n t > for both positive and negative e. 



±(A(2 
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Following j.JSS , let us next define the Priifer variables with some care. For EgR, let u E 
be the sequence of real numbers given via the recurrence relation H w u E = Eu E and the initial 
condition u E (— 1) = sin(#°) and t w (0)w £ '(0) = cos(6 )0 ) where t u (n) = —(n\H u \n — 1). The free 
Priifer phases 8® ,E (n) and amplitudes R^ E (n) > are now defined by 

pO,-E / \ ( cos(6°j E (n)) \ / Un)u E (n) \ 

K « [n) \ sm(9°j E (n)) J ~ V u E (n-l) J ' [6) 
the above initial conditions as well as 

~<c>+i)-ew<|- 

The interest will be on energies E = E v + e in the vicinity of the band edge E u , namely |e| < 
e . Associated to the basis change (J2J), the M e -modified Priifer variables (ii^(n), 8^(n)) G R + xM 
will be introduced next. Define a smooth function m e : R — > R with m e (# + 7r) = m e (8) + 7r and 
< Ci < m' e < C 2 < oo, by 



Then set 0£(n) = m e (0°' B " +e (n)) and 

/ cos(^(n)) \ / Un)u E ^(n) \ 

K " [n) \ sin(^(n)) y ~ e V j ' W 

where the dependence on the initial phase 8^(0) = m e (8°) is suppressed. The oscillation theorem 
as proven in jJSSj implies that the IDS close to the band edge E v is given by 

M(E u + e) = - lim -E(fl») , (5) 

7T n^oo n 

the expectation being taken w.r.t. P. If p = S a , this formula gives the IDS M a of the L-periodic 
operator H a . A similar formula allows to express the Lyapunov exponent in terms of the the 
Priifer variables jJSSj . but this will not be used here. 

The M e -modified phase shift dynamics S eyCT {8) (with energy variation e relative to the band 
edge Ey) is defined via the M e -modified Priifer phase with initial condition 8° = 8 by S ejfT (8) = 
8^(L) — L it J\f a (E u ) where ui = ((cr n )nez, I — 0) and o\ = a. Note that it verifies 

^(S)^'e'«,(:«:») t p>o. (a. 

One then obtains a discrete time random dynamical system S™ u on R defined iteratively by: 

S^ifi) = S^ m {S^\8)) , S°J6) = 8 . 
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Figure 1: Schematic plots of the phase shift dynamics S EtU projected on RP(1) = [— |, ~) and 
for A = — 1 in the elliptic (n t > 0), parabolic (e = and n e = 0) and hyperbolic {n e < 0) regime. 



Replacing the Priifer phases in (j5J) by the phase shifts relative to the band edge E u , the IDS in 
the vicinity of E v is given by 



8M{±e) = j - lira i E($£ (0)) , 
where the sign is chosen such that for positive e one enters the spectrum of H u 



(7) 



Comparing (J2J) and (JUJ), one notes that the phase shift dynamics S eyV can be immediately 
calculated: 



-n e cos(6) + A sin(#) 



A(l - K e ) cos(fl) + sin(0) ' 

Alternatively the cotangent may be used. For A = — 1 the curves are plotted in Fig. 1 for 
three different values of K e . Iterating S eyV gives us a discrete time dynamical system on M, 
which due to the periodicity relation may also be regarded as the lift of a dynamical system on 
MP(1) — [— f , |). If K e > 0, there is no fixed point and therefore the dynamics S e ^ is conjugate 
to a rotation. The rotation number can be calculated explicitely, but it is roughly equal to 7r 
over the number of iterations needed to go through one period. For n e < 0, there are two fixed 
points per period, one unstable and another one stable and (globally) attractive so that the 
rotation number is 0. In the parabolic case K t = there is only one fixed point, instable to 
one side and stable to the other, and the rotation number is still 0. What was just described is 
locally simply a saddle node bifurcation. 



3 Proof 

Fig. 2 shows the elliptic dynamics S t)V as well as a second dynamics which is hyperbolic. The 
latter should be thought of as representing those of the hyperbolic iS £)tT , a ^ z/, which is closest 
to the elliptic case. Let us first argue that the hypothesis imply that Fig. 2 is qualitatively 
correct. Indeed, E u is supposed to be a band edge only of H v so that iS £it7 is hyperbolic for all 
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Figure 2: Plot of the dynamics S €jU just in the elliptic regime, and the next closest hyperbolic 
dynamics (dotted curve). This is the relevant sitution for the study of the Lifshitz tails. 

a v as long as |e| < e for some adequately chosen 6q. Furthermore I n ] is an eigenvector 



t 

of MT®" 'M~ x which implies that the fixed point of the parabolic map S^ v is 9 = 0. For e > 0, 
the map <S £ji , is given by shifting the graph of <S ,v m t° the elliptic regime. By hypothesis, 

M _1 ^ q j is not an eigenvector of T^ v for any a ^ u, hence the fixed points of iS ej(T , for a ^ v 

and |e| < eo, are bounded away from 9 = 0. These facts are resumed in Fig. 2. 

Next let us briefly present the main argument qualitatively. According to (JJJ), the IDS is given 
by the mean rotation number, the average being taken w.r.t. the probability measure choosing 
the upper and lower graph in Fig. 2 randomly. Very close to the band edge, the dynamics S eM 
is only slightly in the elliptic regime and many iterations are necessary in order to complete 
one rotation. During most of these iterations, the angle is in a small interval I = [—a, a] close 
to the origin. If at any of these iterations any of the other dynamics is chosen, the angle is 
immediately again outside and to the left of I (this will be the definition of /). Hence the only 
way to go through I and hence complete a rotation is to always choose the dynamics S €jU until 
the angle is to the right of I. This happens with a very small probability which leads to the 
precise form of the Lifshitz tails. 

The proof of the lower bound now goes as follows. The rotation number of S CjU is equal to 

nL5Af u (e). The number of iterations needed to complete one rotation is K = LSJ ^ ^ + 1 

(here [b] denotes the integer part of b G R). Then 5£,(0) > 9 + vr, but S^ 2 (6) < 6 + tt for all 
9 G R. For sake of concreteness, suppose that E v is a lower band edge so that the sign in (J7J) is 
+ . Next let us set m = KN in ((7j) and decompose S^(9) into a telescopic sum: 



111 1 N 

5A/ " (6) = n L K & N E E ^ {9) ~ S ^' 1] ^ • (8) 



n=l 
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With probability p^ one has Oj = v for j = K(n — 1) + 1, . . . , Kn. In this case, one rota- 
tion is completed and hence each summand can be bounded from below by np^ . Elementary 
inequalities now imply the lower bound in (Q). 

In order to prove the upper bound, set a = max{9 < 1 \S^ e (9) < — 9 V |e| < e , a / u} 
and M = inf{m > 1 | S™ v (a) > vr — a V < e < eo}. By construction, the only way to cross 

- M times. This happens with probability 



[— a, a] is to chose a = v at least K 



LSN v {e 



. If this event occurs, the accumulated phase shift is of order 2a which can simply be bounded 
above by ir. Hence using the same decomposition as in (jHJ), but with the different K, one gets 
the bound 5J\f(e) < jj^Pu ■ As M is finite, this implies the upper bound in (JTJ). 



4 Outlook 

The upper bound in the above argument exploits the following fact: the only way to cross the 
critical region / = [—a, a] is by sucessively choosing the favorable branch S e , v . As this is the 
only way, it seems adequate to speak of a deterministic estimate. A more complete analysis of 
the mean rotation number would also have to take into account what happens in the remainder 
S' 1 \/. Obviously that heavily depends on the precise form of the maps S e ^ a as well as the 
probability measure p. 

A situation in which an analysis becomes feasable is perturbation theory. Suppose that 
dynamics S\^ )(T depend on a supplementary small parameter A giving the order of the L°°- 
distance between all the maps <$a,<e,o-- Then the random dynamics in S l \I can be analysed 
perturbatively in A allowing to calucalate the Lipshitz constants perturbatively. This situation 
arises for example in the Anderson model in the weak coupling limit, namely H w = Hq + W w 
with some L-periodic backround operator H and a random potential V u . Then the gaps of 
Hq remain open for sufficiently small A and the IDS in the gap is still given by the gap label 
of Hq. Within this framework it is also possible to study the IDS away from the Lifshitz 
tails. If all the dynamics at a given energy E are elliptic, then the IDS is simply given by 
N(E) = Jdp(a)Af a (E) + 0(X 2 ) (this can be proven along the lines of Sec. 4.5 of jlEE])- The 
situation becomes more interesting at energies where there are both elliptic and hypobolic 
phase shift dynamics. Then the mean rotation number can be calculated perturbatively via 
a classical ruin problem associated to the passage through an (appropriately chosen) interval 
/. The Lifshitz tails are recovered on one extreme (where an error leads to immediate ruin), 
but this allows moreover to control a cascade of large deviation regimes between the Lifshitz 
tails and the band center. By the same techniques the Lyapunov exponent can be computed 
perturbatively, completing thus the results of |PF| Thm. 14.6] and |JSS1 Sec. 4.6]. The work 
giving a detailed analysis corresponding to these ideas is under preparation. 
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